Search with EGRET for a Gamma Ray Line from the Galactic Center 
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We search data from the Energetic Gamma Ray Experiment Telescope (EGRET) for a gamma-ray 
line in the energy range 0.1-10 GeV from the 10° x 10° region around the Galactic center. Our null 
results lead to upper limits to the line flux from the Galactic center. Such lines may have appeared 
if the dark matter in the Galactic halo is composed of weakly-interacting massive particles (WIMPs) 
in the mass range 0.1-10 GeV. For a given dark-matter-halo model, our null search translates to 
upper limits to the WIMP two-photon annihilation cross section as a function of WIMP mass. We 
show that for a toy model in which Majorana WIMPs in this mass range annihilate only to electron- 
positron pairs, these upper limits supersede those derived from measurements of the 511-keV line 
and continuum photons from internal bremsstrahlung at the Galactic center. 

PACS numbers: 95.85.Pw, 95.35.+d, 98.35.Gi 



I. INTRODUCTION 

Weakly-interacting massive particles (WIMPs) provide 
promising candidates for the dark matter in Galactic: ha- 
los [2,12, 3]. The most deeply explored WIMP candidate 
is the neutralino, the lightest superpartner in many su- 
persymmetric extensions of the standard model Al- 
though the favored mass range for neutralinos is usually 
> 10 GeV, there are other WIMP candidates with masses 
in the 0.1-10 GeV range. For example, neutralinos with 
masses as low as 6 GeV are plausible if gaugino unifi- 
cation is not assumed @. Neutralinos with masses as 
low as 100 MeV are plausible in the next-to-minimal su- 
persymmetric standard model (NMSSM) [(| Q. Also, 
scalar and spin- 1/2 particles with masses in the MeV 
range have been considered [|| to explain the 511-keV 
gamma-ray line observed by INTEGRAL 0, [l(|, a line 
whose strength, as explained in Ref. Q, has defied easy 
explanation from traditional astrophysics. 

One way to detect WIMPs is to search for monoen- 
ergetic gamma rays produced by pair annihilation in the 
Galactic halo [ll| . These gamma rays have energies equal 
to the WIMP mass m x . Such a line spectrum could be 
easily distinguished from the continuum spectrum from 
more prosaic gamma-ray sources (e.g., cosmic-ray spalla- 
tion), and thus serve as a "smoking gun" for dark-matter 
annihilation. 

Since the dark-matter density is highest at the Galac- 
tic center, the flux of WIMP-annihilation photons should 
be greatest from that direction. On the other hand, the 
continuum background should also be highest from the 
Galactic center. We estimate that for a Navarro-Frenk- 
White profile [p|, the WIMP-annihilation flux from the 
10° x 10° region from the Galactic center should ex- 
ceed that from the Galactic anticenter by a factor ~100, 
while the flux of cosmic-ray-induced photons at energies 
O(GeV) is only about 8 times higher from the Galac- 
tic center than from the Galactic anticenter. Thus, the 
Galactic center is the preferred place to look for a WIMP- 
annihilation signal. It is also the location of the 511-keV 
anomaly that has motivated the consideration of lower- 



mass WIMPs. 

In this paper, we search data from the Energetic 
Gamma Ray Experiment Telescope (EGRET) [13| on 
the Compton Gamma Ray Observatory (CGRO) for a 
gamma-ray line in the energy range 100 MeV to 10 GeV 
from a 10° x 10° region around the Galactic center. We 
found no evidence for a gamma-ray line from the Galac- 
tic center in this energy range. From these null results, 



we can bound the cross section 



(TV) 



for WIMP anni- 



hilation to two photons for WIMPs in this mass range. 

The plan of our paper is as follows: In Section HH we 
discuss how EGRET data are cataloged. In Section HTT1 
we reconstruct from the EGRET data the differential flux 
of photons as a function of energy. In Section ITVl we fit 
to the data a model of the flux produced by cosmic rays 
and point sources near the Galactic center. In Section fVl 
we search for a line excess of photons from WIMP annihi- 
lation. In Section fVTl we report upper limits to 
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as 



a function of ra x for WIMPs within the mass range of 0.1 
GeV to 10 GeV for a variety of dark-matter-halo models. 
In Section IVlIl we show that in a toy model in which the 
WIMP annihilates only to electron-positron pairs, this 
upper limit is stronger over this mass range than lim- 
its derived from the 511-keV line and from lower-energy 
continuum gamma rays from internal bremsstrahlung. 



II. SOURCE OF DATA 

We obtained publicly available data from the CGRO 
Science Support Center (COSSC). 1 We used the EGRET 
photon lists (QVP files), which contain event lists of all 
photons detected during a given viewing period. The 
data that we used from these files are the photon's Galac- 
tic latitude, Galactic longitude, zenith angle, energy, and 
energy uncertainty. We also required the exposure files, 
which contain the detector's effective area multiplied by 



1 http: / / cossc.gsfc.nasa.gov/docs / cgro / cossc / egret / 
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TABLE I: Viewing periods used in analysis. The more dom- 
inant viewing periods are in bold and have an exposure of 
> 10 6 cm 2 s sr at 150-300 MeV, over our region of interest. 
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the viewing time of the detector for a particular view- 
ing period multiplied by EGRET's 1-sr held of view. 
The exposure is provided as a function of latitude, lon- 
gitude, and energy range. We also obtained the counts 
files, which contain the number of photons at various 
spatial coordinates and energy ranges within a viewing 
period. The energy bins, along with their respective en- 
ergy ranges, are shown on the COSSC site. 



III. CONSTRUCTION OF GAMMA RAY FLUX 

We begin by constructing the photon differential flux 
as a function of energy. We use data only from a square 
region on the sky from —5° to 5° Galactic longitude and 
—5° to 5° Galactic latitude. Each viewing period covers 
a particular region of the sky, and there were 34 viewing 
periods for our region of interest. These viewing periods 
were found using Table 1 in the Third EGRET Catalog 
PH and are listed in Table [fl 

The differential photon flux can be determined from 
the counts files provided by EGRET, but these provide 
only counts in 10 energy bins, each with a width com- 
parable to the photon energy in that bin. However, we 
will below search for lines with energies spanning the full 
energy range. This analysis is performed (as discussed 
below) by fitting the measured photon distribution to a 
continuum plus a line broadened by a Gaussian, consis- 
tent with the instrumental resolution, about each central 
line energy. We therefore work with the EGRET events 
and exposure files, which list an energy and effective ex- 
posure, respectively, for each photon, and reconstruct the 
differential energy flux in 119 energy bins. Before do- 
ing so, however, we first construct the differential energy 
flux from the events files with the same 10 bins as in the 
EGRET counts files, to be sure that our event-file analy- 
sis recovers the EGRET counts files, the most commonly 
used EGRET data product. 

We first split the data into the 10 energy bins used by 
EGRET. Since the exposure files record a photon index 
value of 2.1 for the photon distribution (a value more-or- 
less consistent with the fluxes arrived at in Figs. [1] and 
[2]), the average energy E avg of photons in an energy bin 
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Variation of the photon index values over the range 
[1.7,2.7] only changes E avg by ~1% for these energy bins 
and by ~0.01% for the 119 smaller energy bins. This 
variation also only changes the average exposures by less 
than 10%, which does not affect our final results signif- 
icantly. Thus, our assumption of a value of 2.1 for the 
photon index is a reasonable one. 

We calculate the differential flux (photons cm 
sr -1 MeV -1 ) from the counts files using 
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e(Ei)AEi 
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where Ei is the average energy of one of the ten large 
energy bins, n(Ei) is the number of photons within that 
energy bin, e(Ei) is the total exposure from the expo- 
sure files over the viewing region within that energy bin, 
and AEi is the size of the energy bin. The quantities 
n(Ei) and s(Ei) are both summed over all viewing peri- 
ods and all positions within the region of interest. The 
uncertainty ap{Ei) in the flux is 
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We assume Gaussian errors in the photon energy. The 
energy uncertainty is just the median of the energy un- 
certainties of the individual photons within that energy 
bin, taken from the events data. 

We then constructed from the events file the photon 
number n(Et) in each counts-file energy bin. We found 
that in order to reproduce the counts data from the 
events file, we needed to reject photons with zenith angles 
greater than 100° and energy uncertainties greater than 
40% of the photon energy. This zenith cut also rejects 
albedo gamma rays from the Earth's atmosphere. The 
photon differential fluxes obtained from both the counts 
files, and the events files (binned in the same way as the 
counts files) are shown in Fig. [T] We were not able to 
match the counts- and events-file photon numbers at the 
first energy bin to within 25%. However, for reasons dis- 
cussed below, we discarded this energy bin (below 0.1 
GeV) from our analysis. 

We then proceeded to construct the differential flux 
from the events files, applying the same photon cuts, 
with narrower bins, to facilitate the analysis in SectionfVl 
We split the data into 119 energy bins, with each bin 
ranging in energy from -E m i n .i = 30 x 1.05 1 MeV to 
£-max,i = 30 x 1.05 l+1 MeV, where i ranges from to 
118. To calculate exposures, we interpolated log[£(.E n )] 
over log(i? n ), where E n is an average energy for a large 
energy bin n, and s(E n ) is the same exposure for the 
large bin n used for the ten large bins earlier. Fig. 14 of 
Ref. [13j shows that the exposures do not vary rapidly for 
energies > 0.1 GeV, and so this interpolation should be 
sufficient for our purposes. The flux is shown in Fig. [2] 
We note that Figs. [1] and agree with EGRET's mea- 
surement of the diffuse gamma-ray spectrum in the same 
region of sky, shown in Fig. 4 of Ref. . We also note a 



3 



r 10- 1 — — — i 

> = = 

<t> : : 

2 r ! A - , <> Counts Files 

en - htftH A Events Files I 

m i m i ; 

E r - 

: , m | : 
to 

1 10 " 8 f 1 

o : : 

Q- r i — — =; 

| 10-° t „^ 

10' 10 2 10 3 10" 

Energy (MeV) 

FIG. 1: The differential flux within ten energy bins with error 
bars denoting energy uncertainty for events data and half-bin 
sizes for counts data. 
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FIG. 2: The photon differential flux using 120 energy bins. 

bump in the differential flux in Fig. [5] at around 3 GeV. 
We believe this artifact is due to the miscalibration of 
Class B photon events (l6j . 



IV. DETERMINATION OF CONTINUUM 
GAMMA-RAY FLUX 

The line we seek is an excess over a continuum, and 
we must therefore model that continuum before we can 
search for an excess. Our aim in this Section is thus to 
find a simple functional form that accurately models the 
continuum over the resolution scales of the instrument. 
A simple linear interpolation over each space of several 
energy-resolution elements would be sufficient, but we 
instead consider several astrophysically motivated func- 
tional forms, although the details of the precise astro- 
physical origin for the continuum are not important for 
our search for a line excess. 

We were able to find a good fit to the continuum by 
a linear combination of three astrophysical sources for 
the diffuse gamma-ray background from the Galaxy. In 
the first source, nuclear interactions, cosmic rays col- 
lide with nuclei in interstellar matter to produce neu- 
tral pions, which decay mostly into gamma rays [l?! ]. 
The second process is bremsstrahlung from cosmic-ray 
electrons interacting with interstellar matter jl7]. The 
third, interior-point-source emission, comes from unre- 
solved point sources within our Galaxy, such as gamma- 
ray pulsars 18] . We also considered exterior-point-source 
emission [19j and inverse- Compton scattering of inter- 
stellar radiation from cosmic-ray electrons, but found 
that the first three sources listed above were sufficient to 
model the flux. Ref. gives the differential gamma-ray 
production functions for the nuclear and bremsstrahlung 
contributions. The production functions are for the 
cosmic-ray spectrum in the solar neighborhood. We as- 
sumed the production functions at the Galactic center 
are proportional to the production functions in the solar 
neighborhood. 

The functional form of the differential flux to which we 
fitted the data was F &t (E) = aF nuc (E) + f3F hrcm (E) + 
aF int (E), where F nuc (E), F hrem (E), and F int (E) are 
the differential photon fluxes from nuclear interactions, 
bremsstrahlung, and interior point sources, respectively, 
and a, (3, and a are amplitudes determined by fitting the 
data. The source functions for nuclear interactions and 
bremsstrahlung are 
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We assume interior point sources to be gamma-ray pul- 



sars. Three pulsars seen by EGRET were the Crab, 
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FIG. 3: The measured and model gamma-ray flux along with 
contributions from nuclear interactions (nuc), bremsstrahlung 
(brem), and interior point sources (int). 



Geminga, and Vela pulsars, which have photon indices 
of -2.12, -1.42, and -1.62, respectively pj]. We ap- 
proximate the photon index as having the average value 
of —1.7, so that 
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The fitted flux [Ffn(Ei)] and the subsequent contribu- 
tions from each physical process are shown in Fig. [3] 



V. ANALYSIS OF EXCESS PHOTONS IN 
GAMMA RAY SPECTRUM 

We next construct a residual number of counts by sub- 
tracting the fitted number iV fit (i?j) = F^ t i K E i )e{E i )/S.E i 
from the observed number N{Ej) of counts. The counts 
N(Ei) and N Rt (Ei) are displayed in Fig.g] 

We take the residual spectrum to be the upper limit 
to the number of photons in each energy bin that could 
come from WIMP annihilation. However, to search for 
the signal we must take into account the finite energy 
resolution. With infinite energy resolution, the WIMP- 
annihilation excess would appear as a monochromatic 
peak over a smooth background distribution. However, 
because of energy uncertainties, each photon captured 
by EGRET will appear to have an energy equal to its 
true energy plus an error, which we take to be Gaussian. 
Thus, monochromatic photons will be spread over neigh- 
boring energy bins. Because our bins are logarithmically 
spaced, the Gaussian will appear skewed, but it will still 
be distinguishable from the background spectrum. 

Suppose our true spectrum before measurement con- 
sists of a continuum C(-Ej) produced by background ra- 
diation and an excess N p of photons with energy E p . Af- 
ter measurement, the continuum will change shape but 



FIG. 4: The spectrum of actual counts, N(Ei), and the fitted 
spectrum, Nat(Ei). 



remain smooth, while the excess will spread out as a 
Gaussian profile over multiple bins. The Gaussian skews 
negligibly, so we approximate the excess as a standard 
Gaussian. Thus, we model the data D(Ei) as 

D{Ek) = C{Ei) + NpfpiEi), (7) 

where f p (Ei) is a normalized Gaussian of the form, 
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In Eq. ([5]), the denominator is summed over all energy 
bins within 3a e p of the Gaussian central energy E p . The 
energy uncertainty <je v at energy E p , is given by 
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where R(E p ) is the dimensionless resolution at energy E p . 
The fractional full width at half-maximum (% FWHM), 
or \/2 In 2 times twice the reciprocal of the resolution, is 
shown for various energies in Fig. 20 in Rcf. [13]. From 
the % FWHM, we produce a table of resolution vs. en- 
ergy, shown in Table [H] We calculate the resolution at 
each energy by interpolating log[R(E)} over \og(E). Be- 
cause the first value for R given in Table |TT] is for energy 
E = 100 MeV, we cannot extrapolate log(i?) to lower en- 
ergies with certainty. Therefore, we restrict our analysis 
to the energy interval 0.1 GeV-10 GeV. 

The number N p (Ei) can be deduced at each energy 
bin in the spectrum by solving Eq. ([7]) for N p , assuming 
D(Ei), C(Ei), and f p (Ei) are known. Each N p (Ei) has 
an uncertainty, 
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TABLE II: Dimensionless resolution R of EGRET at various 
energies. 
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FIG. 5: Ratio of excess photons to the excess uncertainty. 

due to continuum fluctuations. Most bins in the spec- 
trum contain large numbers of photons. Therefore, we 
average N p using gaussian statistics to calculate N p and 
rjjY ■ the value and uncertainty of the excess, for each 
energy bin E p greater than 100 MeV. The resulting ratio 
of N p to a-^ is shown in Fig. O 

Fig. [5] does show statistically significant deviations of 
the data from our model for the continuum. To determine 
if this residual favors the Gaussian model, we compare \ 2 
for a Gaussian model to x 2 for a constant-excess model. 
We calculate \ 2 f° r both models over a ±3cte range 
centered at the excess center. The Gaussian is N p f p . 
We also compare the residual with a constant excess iV c , 
where dN c /dE is constant and N c is proportional to the 
energy-bin size. We normalize N c such that the lowest 
energy bin ?><Je v from the Gaussian center has 10 pho- 
tons. We compared \ 2 f° r the excess at energies E = 210 
MeV and E = 2000 MeV, two energies that have high ex- 
cess photons to excess uncertainty ratios (see Fig. [3]). At 
both energies we found % 2 to be smaller for the constant 
excess, a simpler model, than for the Gaussian. Thus, 
we show that the residual does not favor the Gaussian 
model, and we do not attribute any of these deviations 
to a WIMP-annihilation line (see Fig.[S]). Rather, it ap- 
pears that there is some continuum contribution that our 
analysis has not taken into account. 

We therefore use N p to calculate an upper limit to the 
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FIG. 6: The residual number of counts (crosses) and the ex- 
pected Gaussian (solid curve) from a smeared line excess. The 
top panel shows the residual and Gaussian at E — 210 MeV, 
while the bottom panel shows the same at E — 2000 MeV. 
The ratio of excess photons to excess uncertainty is high at 
these energies. Notice in both panels the residual does not 
resemble the Gaussian. For E = 210 MeV and E = 2000 
MeV, respectively, \ 2 f° r the Gaussian is 17.3 and 38.0 and 
X 2 for the constant excess is 11.0 and 23.9. Thus, the Gaus- 
sian model is not favored. 

line flux. This line flux is different from the differential 
flux used in previous sections in that this flux is not di- 
vided by the energy bin size. Since N p has positive and 
negative values, we take the 2a upper limit to the line 
flux ^ u (E p ) to be 

The 2<T upper limit to the line flux is shown in Fig. [71 

We illustrate the reliability of the upper limit to the 
line flux by repeating the analysis in Section [V] for a 
sliding-window continuum model. At each energy bin 
Ei we fitted the diffuse flux data within 3 to 9er E . of 
Ei to a single power law. The amplitude and index of 
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TABLE III: Parameters for each profile type. 
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FIG. 7: Upper limits to the line flux $ u from the Galactic 
center. The solid line is the upper limit derived from the 
continuum model in Section llVl The dashed line is the upper 
limit derived from the sliding window technique. 



the power law, which varied with energy bin, were then 
used to construct the background radiation continuum 
C(Ei) in Section [V] needed to search for a line excess. No 
significant excess was found, and an upper limit to the 
line flux was determined. This 2a upper limit, shown in 
Fig. agrees quite well with the previous upper limit 
in Section [V] except around 3 GeV, where the previous 
upper bound is more conservative. To be conservative, 
we chose the upper limit to the line flux from the multi- 
component continuum fits, for the rest of our analysis. 



VI. UPPER LIMITS TO THE ANNIHILATION 
CROSS SECTION 

If WIMPs comprise the Galactic halo, then the flux of 
line photons from WIMP annihilation is (for Majorana 
WIMPs) 



{(TV 



p\dl, 



(12) 



where $ is the line flux of photons in units of photons 
cm~ 2 s _1 sr _1 , (eru) 77 is the velocity-averaged cross sec- 
tion for the WIMP to annihilate to two photons, m x is 
the WIMP mass (which is equal to the photon energy 
E 7 ), and p x is the density profile of the WIMP halo. 
The integral is along the line of sight, and dl is the dif- 
ferential distance along the line of sight. The residual in 
the previous Section gives the average line-of-sight line 
flux within a 10° x 10° region around the Galactic cen- 
ter. Therefore, we integrate Eq. (|12j) over our viewing 
region to find the relation between (ert;) 77 and m x . 

The density profile of the WIMP halo must be known 
in order to integrate Eq. (JT3J) . The functional form of 
the halo density profile is motivated by theory and simu- 
lations, with parameters chosen for consistency with the 
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FIG. 8: The 2a upper limits to the velocity-averaged annihi- 
lation cross section (au) 77 as a function of WIMP mass for 
various halo-density profiles. 



measured Milky Way rotation curve. We assume the fol- 
lowing parametrization of the density profile, 



P(r) = Pa 



(ro/aF[l + (ro/a) a ]V-iV a 
{r/ap[l + (r/a) a ]<P-i)/ a ' 



(13) 



Here, po is the local density of the halo at the Solar Sys- 
tem; ro is the distance from the Solar System to the 
Galactic center, which we take to be 8.5 kpc; a is the 
core radius; and a, /?, and 7 are parameters that deter- 
mine the halo model. Various combinations of a, 0, and 
7 have been used in simulations and are of particular 
interest. We chose to study the Ka and Kb profiles pro- 
posed by Kravtsov et al. [2(|; the NFW profile proposed 
by Navarro, Frenk, and White and the modified 

isothermal profile, or Iso, which is commonly used. These 
profiles are listed in Table IIIII The quantities po and a 
are chosen for each profile so that the profile will account 
for the Galactic rotation curve. These values are taken 
from Fig. 5 in Ref. [ll|. We insert each of these profiles 
into Eq. (|12|) and integrate over our viewing region to 
find the line flux $ in terms of (av)jj and m x . 

The resulting upper limit to the annihilation cross sec- 
tion (cru) 77 is shown in Fig. [5] as a function of WIMP 
mass m x for each halo model listed in Table [TTTT 
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VII. DISCUSSION 

To illustrate the possible utility of this new bound, 
we consider a toy model in which WIMPs are Majorana 
fermions that couple to electrons via exchange of a scalar 
boson (the U boson dHHU) of mass m\j (assumed to be 
much heavier than both WIMPs and electrons) through 
the Lagrangian density, 
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where CVy and /a e are axial couplings of the U boson to 
the WIMP field \ and the electron field ip e , respectively. 
Annihilation of WIMPs with O(MeV) masses to electron- 
positron pairs has been considered as a possible explana- 
tion § for the observed flux, <& 5 n = 9.9^2;! x 10~ 4 pho- 
tons cm -2 s _1 @, of 5 11-keV photons as measured at the 
Galactic center by the SPI camera on the INTEGRAL 
satellite. In this scenario, positrons from WIMP annihi- 
lation then annihilate with electrons in the IGM to pro- 
duce these 511-keV photons. The annihilation rate — and 
therefore the cross section for annihilation to electron- 
positron pairs and thus the coupling Gu fAe/mfj — are 
determined by the flux of 511-keV photons. More pre- 
cisely, the 511-keV flux determines an upper bound to 
this annihilation rate, cross section, and coupling, but 
we will here suppose the entire 511-keV flux to be from 
positrons from WIMP annihilation. 

Ref. 124[ pointed out that if WIMPs annihilate to 



electron-positron pairs, they can also undergo annihi- 
lation to an electron-positron-photon three-body final 
state, a process we refer to as internal bremsstrahlung. If 
(uv) e + e ~ is the cross section for annihilation to electron- 
positron pairs (as calculated, e.g., in Refs. [22l I23I |25| ) . 
then the differential cross section for bremsstrahlung of 
a photon of energy _E 7 is 
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4m x (m x — -E 7 ), and a ( 



(15) 
is the 



where s 

fine-structure constant. The quantity E^d(av)Br/dE~ 1 in- 
creases roughly linearly with E-y for E 7 < m x and peaks 
at a value (for our WIMP mass range of 0.1-10 GeV) less 
than 10% smaller than the WIMP mass. The measured 
upper limits to the flux were approximated in Ref. pij 
E^Bx/dE^ < 7 x 10~ 3 MeV cm" 2 s" 1 sr" 1 over the 
energy range 1-100 MeV. This flux was averaged over a 
region on the sky centered at the Galactic center from 
-30° to 30° Galactic longitude and -5° to 5° Galactic 
latitude. For the purposes of this illustrative exercise, we 
extend this bound up to 10 GeV (roughly consistent with 
the line limit we have derived). 

Each annihilation to an electron-positron pair pro- 
duces two 511-keV photons either directly (7% of all an- 
nihilations) or by producing positronium and decaying 
(23.3% of all annihilations); the rest produce noncon- 



tributing continuum photons [24|, |26j . The resulting flux 
of 511-keV photons is (for Majorana particles) 



$511 = 



47rm 2 



p 2 Y dl <m, 



(16) 



where £ = 0.303 is the fraction of positrons that undergo 
two-photon annihilation, the dl integral is along the line 
of sight and the dfl integral is over the SPI camera's 
field of view, a 16°-diameter circle around the Galactic 
center. Likewise, the differential flux of photons from 
internal bremsstrahlung is 



dEr, 



d(av)Br/dE 1 
87rm 2 Aft 



p 2 dl dSl, 



(17) 



where Aft ~ 0.182 sr is the solid angle over the 60° by 
10° Galactic region mentioned earlier. 

The two-photon annihilation cross section (a"i;) 77 for 
the Lagrangian of Eq. f|14|) is given by [27j 



«^fie |/fe)|2) 



S 4 

ir^mjj 



(18) 



where £ e = to 2 /to 2 , J(f e ) = \ [1 + £ e J(f e )], and J(£ e ) is 
given by 



.1 1 + yi - £ e in 



(19) 



for £ e < 1. For our WIMP mass range 0.1-10 GeV, £ e < 
1 and I(£ e ) — 1/2. The cross section for annihilation to 
electron-positron pairs (av) e + e - is given by [23l. [25| 



{(TV 



2nmfj 



4 2- 

-m x v\ 



(20) 



where u 2 = |t> 2 is the mean-square center-of-mass ve- 
locity and v c ~ 220 km/s is the WIMP rotation speed, 
assuming the electron energy E e = m x 3> m e and 
mu ^> m x . We use Eqs. IT8|) and (|20l to derive up- 
per limits to the coupling Cu fAe/mjj appearing in the 
Lagrangian of Eq. (TTJJ . 

Fig. [9] shows the upper limit, assuming an NFW 
halo-density profile, to the coupling CufAe/mfj from 
measurements of the 511-keV line the limit to the 
bremsstrahlung-photon flux [2~ij , and our la limit to the 
line-photon flux. We see that for the model assumptions 
and WIMP mass range considered here, the limit to the 
two-photon annihilation cross section derived from our 
2cr limit to the line-photon flux is the strongest of these 
three. At first, this result may seem surprising, given that 
the two-photon annihilation process is higher order in a e , 
but this suppression is counteracted by the helicity sup- 
pression of the cross section for annihilation of Majorana 
fermions to electron-positron pairs. Refs. [H, [2£| con- 
sidered also gamma-rays from in-flight annihilation from 
c + c~ pairs, but their analysis was restricted to energies 
< 100 MeV. 
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FIG. 9: Upper limits to the ratio Cu fAe/mfj as a function 
of WIMP mass for the NFW halo-density profile. The limits 
were calculated from the observed 511 keV emission, the con- 
straints on internal bremsstrahlung, and our derived limit to 
the line photon flux. 



Of course, the 2a limit to the line-photon flux may 
not always provide the best limit to the two-photon an- 
nihilation cross section for every WIMP model. It may 
well be that other models — e.g. , those in which the dark- 
matter particle is a scalar [30| — can produce a ratio of 
511-keV photons to line photons large enough to cause 
the 511-keV limit to supersede the line photon limit. 
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